Motion by Volume Preserving Mean Curvature 
Flow Near Cylinders 

David Hartley* 
Abstract 

Center manifold analysis can be used in order to investigate the stability 
of the stationary solutions of various PDEs. This can be done by considering 
the PDE as an ODE between certain Banach spaces and linearising about the 
stationary solution. Here we investigate the volume preserving mean curva- 
ture flow using such a technique. We will consider surfaces with boundary 
contained within two parallel planes such that the surface meets these planes 
orthogonally. The stationary solution we consider is a cylinder. We will find 
that for initial surfaces that are sufficiently close to a cylinder the flow will 
exist for all time and converge to a cylinder exponentially. In particular, we 
show that there exists global solutions to the flow that converge to a cylinder, 
which are initially non-axially symmetric. A similar case where the initial 
surfaces are compact without boundary has previously been investigated by 
Escher and Simonett |0. 

1 Introduction 

Given a sufficiently smooth, compact initial hypersurface Oo = Xo(M n ) c R' ,+1 we 
are interested in finding a family of embeddings X : M" x [0, T) — > R w+1 such that 

dX 1 r 

— = (h-H)vn„ X(-,0) = Z , h = — Hdfi t , (1.1) 

or \ll t \ J M n 

where H = YI! = \ K i i s tne mean curvature and /q are the principal curvatures of the 
surface Q. t = X(M n , t) = X t (M n ), vn, and dfi t are a choice of unit normal and the 
induced measure of £l t respectively. Note that the flow (ll.lt decreases surface area 
while preserving volume (see Section 01. 

For convex hypersurfaces without boundary it is known that a unique solution 
exists for all time and that the hypersurfaces converge to a sphere as t —> oo Q. 
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Other results include average mean convex surfaces with initially small traceless 
second fundamental form converging to spheres Q and surfaces that are graphs 
over spheres with a height function close to zero, in a certain function space, con- 
verging to spheres The case where the initial surface has a boundary has also 
been studied. In this case it is assumed that Q is smoothly embedded in the domain 

W = [x e R M+1 : < x n+ i < d) , 

with d > and d£lo cz dW. The volume enclosed by Q.q and W will be labelled V 
and vn f is chosen to point outward. The boundary conditions for the flow are then 
that Qo intersects dW orthogonally. Assuming Qo to be rotationally symmetric it 
was proved in [Q] that the flow exists for all time and converges to a cylinder in W 
of volume V under the assumption 

Ifiol < V -i- (1-2) 
a 

This constraint is needed to ensure that the solution never touches the axis of ro- 
tation, so that no singularities develop. Short time existence of the flow can be 
proven using standard parabolic techniques for initially smooth surfaces, however 
in this paper we are concerned with a larger class of initial surfaces so need a more 
general existence theorem (see Section |2). 
The main result of this paper is 

Theorem 1.1. If Qo is a graph over the cylinder ^Si, such that R > d ^~^ and 
with height function sufficiently small in h ^ \^r\ < /? < 1, then its flow by 

dz 

di.il ) exists for all time and converges exponentially fast to a cylinder as t — > oo, 
with respect to the C k topology for any k e N. 

An important point to note is that the limiting cylinder may not be the same as 
the cylinder ^ , thus the cylinders are stable critical points of the PDE (11.11 ). but 
not necessarily asymptotically stable. An immediate consequence of this theorem 
is that there are non-axially symmetric surfaces that converge to a cylinder under 
the flow (ll.ll ). since there are non-axially symmetric surfaces in any h x ^ a neigh- 

Tz 

bourhood of a cylinder. Lastly note that the condition R > i^JEEL j s Q f a similar 
nature to the condition (11.21 ); in fact for a cylinder dl.21 reduces to R > nd. In view 
of this, when considering the flow of axially-symmetric surfaces that are close to a 
cylinder we reproduce the results of |[T] under a weaker condition. The condition 
is also sharp, in the sense that if R = i^El then there is a one-parameter family of 
unduloids, satisfying the boundary conditions, such that ff" , is the only cylinder in 
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the family, i.e., there are stationary solutions to the flow arbitarily close to , that 
are not cylinders. 

In Section|2]of this paper we convert the flow (11.11 ) to a PDE for the graph func- 
tion and also introduce the function spaces and notation that will be used through- 
out the paper. The section ends with some important theorems about the flow. In 
Section [3] we consider the problem as an ODE on Banach spaces and determine 
the spectrum of the linearisation of the speed. This then allows us to find a center 
manifold of the flow, which is seen to consist locally entirely of cylinders. Lastly 
exponential attractivity to this manifold is proven. 

The author would like to thank Maria Athanassenas for her initial suggestion of 
working with cylinders, support, advice and help in preparing this paper, and also 
Todd Oliynyk for his advice and encouragement in dealing with center manifolds. 
Finally to Monash University and the School of Mathematical Sciences for their 
support. 



2 Notation and Preliminaries 



In this paper we consider M" 



Xo(p) = P + Po(p)v^ d (p), p € so that X is a graph over ffg d . 



^ a given cylinder of radius R and length d, and 

The volume 

form on such a graph surface will be denoted by dp p and we let p p be the function 
such that dfip = p p dpo. We use the notation q e ^S -1 , the (n - l)-dimensional 
sphere of radius R, and z £ [0, d], so that (q,z) £ ^ra- Proceeding as Escher 
and Simonett |3] and converting the flow to an evolution equation for the height 
function p : ff" , X [0, T) — > R it is known that up to a tangential diffeomorphism 
the flow (11.11 ) is equivalent to 



=0,d 



= 0, pG,0)=p , 



(2.1) 



where H p = H p (-, t) is the mean curvature of the surface defined by p(-, t) and g p 
denotes the inverse metric on ^ 1 +p d (see JH). 

The graph functions p are chosen in the little-Holder spaces, which are defined 
for a € (0, 1), k € N, an open set U c W and a multi-index f3 = (fit,... ,j3„) with 
10 = £" =1 jS; as follows (see D): 



h a (D) = \pe C a (V) : lim sup 

[ x,yeU, 0<\x-y\<r 

h k+a (U) = {pe C k+a (0) :lfpeh a (u) for all J3, \/3\ = k) , 



\p(x)-p(y)\ 

\*-y\ a 
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where D is the derivative operator on R" and C a , C k+a are the Holder spaces: 

\p(x)- P (y)\ 



C a (D) = lpeC(u): sup 

{ x,yeU,xJ=y 



< oo 

\x-y\ a 



\\p\\h^(u) = \\pWc^{u) = Mlc*(t?) + 2 s " p | v V |.. 



C k+a (jj) = {peC k (U) :LfpeC a (U) for all (3, ]0\ = k] . 
The norms on the little-Holder spaces are the same as those on the Holder spaces, 

\DPp{x) - DPp{y)\ 

Tfc x,yeU,xj=y 

where 

\\P\\ C *(U) = 2 SU P KpW| ■ 

\P\<k 

The little-Holder spaces can be extended to a manifold by means of the atlas of ^ ,. 
In addition, it is known that the little-Holder spaces are the continuous interpolation 
spaces between the C k spaces (see (H), 

h ke (^ d ) = {c(^ d )x k {^ d )) d , 

where (•, •) is an interpolation functor for each 9 € (0, 1) and is defined for Y c X 
as follows: 

(X, Y) B = {xeX: lim t~ e K(t, x, X, Y) = ol , K(t, x, X, Y) = inf (\\x - a\\ x + t \\a\\ ¥ ) . 

The interpolation space h a ( ^ can also be viewed as the closure of the smooth 
functions on the cylinder, C°° (j^g d J, inside the space of Holder continuous func- 
tions C"^^^ (see [8 ]). To simplify the notation involving the boundary condi- 
tions, let Y ytfgd) represent any function space on the cylinder and set 



F-(^) = {peF(^):| 



-0 

z=0,d 



dz 

With this notation we also have that 

/ at 2+a( n \\ i^^Rd) 9 ^ ~T 

\h a (% d ) , h r (^id) i = I h a i 2e fe) 6>±f-,e± i-f ' (2 - 2) 

where a e (0, 1) - a result that follows from the Reiteration Theorem ([8 ] Theorem 
1.2.15) along with the characterisation of interpolation spaces in Chapter 3 of JH. 

For an operator between function spaces G : Y — > Y we denote the Frechet 
derivative by dG. A linear operator, A : Y c X — > X, is called sectorial if there 
exist 6 6 (f ,7rJ and weR and M > such that 
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(i) p(A) D = [A € C : ^ * to, || arg(i - oj)\ < 9], 

(ii) ||i?(i, A)\\ L{X) < ^ f or a n X e 5 e , w , 

here p(A) is the resolvent set, i?(/l,A) = (AI - A) -1 is the resolvent operator and 
II ' IIl(Z) is the standard linear operator norm (see |)8]). 
The following are some facts about the flow. 

Lemma 2.1. For an initially smooth, compact hypersurface with boundary con- 
tained in dW, the flow di.il) together with Neumann boundary condition: 

( i) preserves the enclosed volume V, 

(ii) decreases the surface area |Q t |. 

Proof. Both parts are covered in |fl]. Using the divergence theorem and setting 
E, c W to be the (n + l)-dimensional space with boundary Q. t we have the first part 
of the lemma as follows 

— = I div — dx = I — • v dEl dS = I (h - H) dp. t = 0, 
at J El ot J dEt dt J M n 

where the integral over the part of the boundary contained in dW vanishes due to 
the boundary conditions. The second statement is through a direct calculation of 
the evolution equation for the area 

—— = - \ {H-hf dp,, 
at Jm" 

□ 

Note that Lemma |2T1 is also true for hypersurfaces without boundary. We now 
give a general short-time existence theorem 

Theorem 2.2. For < [3 < 1 and a given hypersurface Clo of class h l+ @, the 
flow di.il ) has a unique classical solution {Q t : t € [0, T)\ for some T > 0. Each 
hypersurface D. t is of class C°° for t € (0, T). Moreover, the mapping t — > Q. t is 
continuous on [0, T) with respect to the h x ^ -topology and smooth on (0, T) with 
respect to the C°° -topology. 

Proof. This result follows from the main result by Escher and Simonett in 0. It is 
clear that the arguments in their proof still hold with h 2+a , h l+ @ and h l+ @° replaced 
by h 2 g a , and h\^° respectively due to equation (12.2b - □ 

3c 3c 3c 
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3 Stability around Cylinders 

The flow in equation (12.11 ) can be considered as an ordinary differential equation 
between Banach spaces. SetO<a</?< 1 and define 



G : h 2 ; a h a (% d ) , G(p) := L p (h p - H p ) , L p := ^1 + gWipVjp. 

The flow (12.11) is then rewritten as 

p'{t) = G(p(t)), 0<t<T, p(0) - p„ € h l ; p (^i d ) . 



(3.1) 
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Lemma 3.1. For the linearisation, dG(0), ofG it holds 



in-l \ n-l 

8G(0)U = I — — + tSfgn I U 



R 2 ®R4 R 2 

forueh 2 ; a (^ d ). 

Proof. Firstly note that L = 1 and that dL p \ p=Q = 0. Hence for u e h 2 ^ a ^ d ) 



dh \ n u = d 



1 



p=0 



r 






I duo d 


\ H p n p duo 


u 


Vn 4 


) 


p=0 







\ \i P duo 


u 




p=Q , 



= i I, f dH p \ u + H Q dn p \ ud[iQ-H Q \ dfi p \ udfi Q 

pj^J W% U6 R,d ) 

= f 9H p \ udno. 



It was shown in [4] that 



dH, 



P\p=0 



so combining these results gives, for u e h 2 f a A 



dG(0)u = I n —±- + A VS \ u -f I — - + A^n )w(i^o- (3.2) 



4, ( R 2 



R 2 *R4 J v „ \ R 2 

The divergence theorem and the boundary conditions give the result. □ 

As we are considering the flow locally about p = it is convenient to rewrite 
(13.11 ) highlighting the dominant linear part 

p\t) = dG(0)p(t) + G {p{t)) , G (u) := G (u) - dG(0)u. (3.3) 

Lemma 3.2. If R > d ^T^ then the spectrum cr (<9G(0)) of 3G(0) consists of a 
sequence of isolated non-positive eigenvalues where the multiplicity of the eigen- 
value is n + 1. 

Proof First note that dG(0) : h 2 + a (^^ -» ^ (^) and since h 2+a (tf£ d ) is 

fiz (iz 

compactly embedded in h a \fo%A the spectrum consists solely of eigenvalues. To 
characterise the spectrum of dG(0) we first look at the spectrum of 

A := — — + Asf* . 

In the following we deal with the natural complexification without distinguishing it 
notationally. The operator A is self adjoint with respect to the L2 inner product on 
h 2 g a y^gj), since the boundary terms disappear for functions in h 2 + a when 

dz ' dz 

integrating by parts. This means that the spectrum consists of real numbers. For an 
eigenfunction v of A assume that v(q,z) - X(q)Z(z), so using A<^ = Ayn-i + |^-: 

1 3 2 \ 

+ Ayn-i + -^\X{q)Z{z) = AX(q)Z(z). 

So, after expanding, 

Z(z)Ayn- 1 X(q) + X(q)Z"(z) + - MX(q)Z(z) = 0. 

Therefore 

Ayn-iX(q) _ z // (z) j n _ l 



X(q) Z(z) \ R 2 

Setting both sides equal to a constant £ e R 

Ayn-iX(q) = -£X(q). 
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The eigenvalues of the Laplacian on are well known to be ~ /(/ ^"~ 2) for I e 

N U {0}, with eigenfunctions the spherical harmonics of order I, Xi p (q) = Y{ p (q), 
1 < p < Mi, where 



Therefore = l ^ 1+ ^ 2 ^ so the other separation equation becomes 

Z"(z) = (a + l ll^±Z^ z(z) , z'(0) = Z' W = 0, 
where we included the boundary conditions to ensure that v e h 2 g +a h?» d ). The 

dz 

eigenvalues and eigenfunctions to the second separation equation are given by 

/mV 1(1 + n- 2) -{n- 1)\ 

+ - ^ -J, m,/eNU(0| 

and 

Z m (z) - cos|^— j. 

So Ai m are the eigenvalues of A with multiplicity M/ and eigenfunctions spanned 
by 

(mnz\ 

vi,p, m (q, z) = cos ^— J 1 < p < Mi. 

Now since the spherical harmonics are dense in the continuous functions on 

and the trigonometric functions are dense in the continuous functions on [0, d], we 

have found every eigenf unction of A. 

Returning to the spectrum of dG(0), vo,i,o = 1 is still an eigenfunction but with 
eigenvalue A = 0. The operator dG(0) is also self adjoint with respect to the L 2 
inner product on h 2 ^ a To see this, consider u, w e h 2 + a (^ d J, then 

dz dz 

X(dG(0)u) w dfjQ = f \Au — T u d/io w dfio 
MA r ) 

= I (Au) w dpo -r u dp.Q I w d^o 

= I uiAw) duo-^—r- I udfiQ-r wd/uo 
M 4 R J K d 

- \ u\Aw-^—^--f wdfio\dfiQ 



X 



"R,d 

u (dG(0)w) dfiQ. 
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Therefore we need only consider other eigenf unctions, orthogonal to vo,i,o = 1, in 
order to characterise the spectrum. This means that for an eigenf unction u 



f 



u d/jQ = 0, 



hence by Lemma 13 - 1 L dG(0)u = Am. The remaining eigenfunctions of 3(7(0) are 
then the remaining eigenfunctions of A, with the same eigenvalues. Therefore for 
m, I e N U {0} and m 2 + I 2 > 1 the eigenvalues are 

(m 2 n 2 l(l + n-2)-(n- 1) 

h,m = ~ — 7, 1- 



d 2 R 2 



with eigenfunctions 



(mnz\ 

vi,p,m(<l> Z) = cos [— J Y i, P (l^ l<P<Mi. 

It is clear from these expressions that = 0, and we relabel the constant eigen- 
function as vi o,o = 1 for convenience, so zero is an eigenvalue of multiplicity 
Mi + 1 = n + 1. If m e N we have A\, m < 0. Also if / > 2 and m e N U then 
Ajj n < 0. The last case is for / = and m > 1 where the eigenvalues are bounded 
above by 

I m 2 n 2 n — 1 \ In 2 n - 1 \ 
\ m = -^—-—j^--—j = A ,. 

These eigenvalues are all negative under the condition R > given in the 

statement of the theorem. □ 

In what follows, we set P to be the projection from h a (^o^^j onto the A = 
eigenspace given by 



Pu :- ^ (w,vi, p ,o)vi,p,o, 



p=0 



where we use (•, •) to denote the L2 inner product on h 2 ^ a Because dG(0) is 

Hz 

self adjoint with respect to this inner product, it is clear that PdG(0)u = dG(0)Pu - 
for every u e h 2 g a Due to this h 2 f a can be split into the subspaces 

dz dz 

X c = P\h 2 l a M and X s = (I - P) lh 2 + a (%%J)\ called the center subspace and 



stable subspace respectively. We are now in a position to apply Theorem 9.2.2 in 
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Theorem 3.3. For any teN, there is a function y e C k 1 (X c , X s ), such that y (k J) 
is Lipschitz continuous, y(0) = <9y(0) = and A\ c = graph(y) is a locally invariant 
manifold of equation \3.3§ with dimension n + 1 provided R > d ^~ l ' . 

Note that by locally invariant it is meant that there exists a neighbourhood 
of zero in A c X c such that if po € graph (y\\) then the solution to (13.31 ) is in 
graph (yU) for all time or until Pp(t) <t A. 

Proof. In order to know that the assumptions of Theorem 9.2.2 in are satisfied 
in our case, we use Equation (I2.2t to show that for any chq satisfying < ao < a < 1 
the following holds 



The assumptions are therefore satisfied by proving there exists a neighbourhood of 
0,Oc /j 2 fl +a J, such that G e C°° (<9, h a J) and, for all p € <9, <9G(p) is 

the part in h a of a sectorial operator A p : /z^ -> h a ° (tfgj). It was 

shown in [2], [4] that the mean curvature can be written in the form 



1 In - 1 



where Q{p) is a second order uniformly elliptic operator that smoothly depends on 
p and its first derivative. This means that G is a smooth operator inside a neigh- 
bourhood, 0\, where if p € Oi we have p(p) > -R for all p e ,. 

As 2(0) = -A%?» rf we have 3// p |^ =0 = -(^r + A<^ rf ), this operator is uni- 
formly elliptic, hence its negative is sectorial. Now the operator 

dz 

In -I \ n-l C 

A u = \- W + A Kd }u- — j-^u dp , 

can be seen to be sectorial: using the definition of sectorial there exists M > 0, 
6 (|, n) and co e R such that 



\A - to\ \\u\\ 2 +aoj x < M 

" d \ J R,dl 



A+ dH p \ p=0 \u 
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for all A e Sg i0) and u e h 2 *" A Therefore 



\U-A Q )u\ 



n- 1 



+ A<^» \\u 



n-l r 



ud/uo 



A + dH, 



P\p=0 



n- 1 



-r u 



> 



n-l. 



n- 1.. 



> \\U ,2+a. 



3z 25 



= w 



1 - T 71-1 



M 



R 2 \A\ j 



|M|l C fe)' 



By the reverse triangle inequality and by taking Re(A) > oj\ := (2|w| + 2M ^" 2 — j we 
get 

IP - Ao)«IU fe) > Ul (1 - ^ - ll»H C o fe) 



7!? 



^ W II II 

- 2^ l|M|l C°fe)' 

dz 

Thus, by applying Proposition 2.1.11 in JSJ, Ao is sectorial. This then implies by 
Proposition 2.4.2 in (H that A p : h 2 ™ (%* d ) -> h a ° (tf£ d ) is sectorial for all p in a 

neighbourhood of zero, 2 c /j 2 fl +Q, ° A Since <9G(p) is the part of A p in /j 0, (<«f" ) 

Bz -\ 

the conditions of Theorem 9.2.2 in [8] are satisfied by setting O = 0\ n 02- n 
We now set 

5:-jp£ /i 2 /" (<^J : graph(p) is a cylinder J . 

Lemma 3.4. M c coincides with the set S on a neighbourhood of zero, provided 
R > 



din- 1 



Proof. Firstly, note that due to Theorem 2.4 in J9), Theorem 2.3 in iTTOl can be 
applied to conclude Af contains all equilibria of (13.31 ) with Ppo e A. Also note, 
the center manifold, Af , is defined differently in iPTOl as compared to [8], however 
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they can be seen to be equal on A. The rest of the proof follows from [3 1, where 
p e S is formulated in terms of eigenf unctions: 



n 

p(y) = ^wi,p,o -R + 
P =i 

with (y\ , . . . , y n , 0) e R M+ 1 being the point of intersection of the axis of rotation with 
the x n+ \ -plane and y$ := R' - R the difference between its radius and the radius of 
. This map is smooth on a neighbourhood U of e K" +1 and its derivative at 
zero is given by 

n 

dp{0)y = j)y*vij,.o. jeIR' ,+1 . 

p=0 

The map taking y to the coordinates of Ppiy) with respect to the basis vi^Oj ^ 
p < n, is then found to have derivative at zero equal to the identity and hence is 
a diffeomorphism from U onto its image, possibly making U smaller. This means 
that the projection of S\u '■= {piy) '■ y £ U] is an open neighbourhood of e X c . 
This can be made to coincide with A (after possible renaming) and since <S|;y c M c 
we conclude that S and M c coincide locally. □ 
We now prove the main result. 

Proof of Theorem 177/1 As in Q we can use the same arguments in (H to prove 
the following result. Given k € N and a> € (0, -Ai), there exists a neighbourhood 
Ui(k, co) of in h 1 ^ with the property: Given po € U\, the solution p(-,po) 

of (13.11 ) exists globally and there exist c(k, oj) > 0, T(k, oj) > 0, and a unique 
Zo(po) £ A such that 

\\p(t,p ) - (zo + 7(zo))\\ C kt V n\ < ce-**\\(I - P)p - vPPoihwfa,,)' 

for t > T; here zo and y(zo) are added in C k ( < ^ (/ )- This proves that p(f) converges 
to an element of Af U, which by the above Lemma is a cylinder. □ 

Corollary 3.5. Let Q.q be a graph over a cylinder, of radius R > d ^ n_1 1 with height 
function po such that the solution, p(t), to the flow f l2.il ) with initial condition po 
exists for all time and converges to zero. Then there exists a neighbourhood, O, of 

po in h l r ^ A < P < 1, such that for every uq € O the solution to d2.il ) with 

55 

initial condition uq exists for all time and converges to a function near zero whose 
graph is a cylinder. 



vp=i 



+ (R + yo) 2 -^]yl 
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Proof. This follows by the same arguments given in J5j. We set U c h 1 ^ (^SJ) to 

be the neighbourhood of zero given in Theorem 11.11 Since p(?) converges to zero 
in the C k topology, there exists a time T such that p(T) € U and as U is open there 
exists an open ball B e {p(T)) c U of radius e centred at p{T). By the finite time 
continuity of the volume preserving mean curvature flow there exists a ball Bg(po) 
such that if uq e Bg(po) then the solution, u(i), to (12.11 ) with initial condition uo 
exists for t e [0, T] and u(T) e B £ {p(T)). Since u(T) is in £/, by Theorem II. li the 
solution to (12.11 ) with initial condition u(T) converges to a function near zero that 
defines a cylinder. By uniqueness of the flow we get the result. □ 
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